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Abstract 

This study investigates the effectiveness of ventilation design strategies for general hospital wards in terms of 
virus removal capacity. A typical semi-enclosed six-bed general ward of Hong Kong hospitals and three 
respiratory viruses, namely Middle East respiratory syndrome coronavirus (MERS-Col/), severe acute respira¬ 
tory syndrome coronavirus (SARS -CoV) and H1N1 influenza virus, were chosen for the computational fluid 
dynamics (CFD) simulation of airflow field and virus dispersion inside the ward. The results demonstrated that 
the location of an infected patient would affect the infection risks to other occupants and healthcare workers 
inside the same hospital ward, and an increased air change rate in the ward could reduce the risk of infection 
from direct contact and inhalation. It was found that an air change rate of 9 h -1 could effectively minimize the 
deposition and floating time of respiratory virus particles while maximizing energy efficiency. This study 
should provide a useful source of reference for the hospital management to mitigate the risk of infection 
with MERS or other airborne transmitted viruses through better ventilation design strategies. 
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Introduction 

According to a statistical report by the Hong Kong 
Hospital Authority, 1 the in-patient discharges and 
deaths were continuously increasing from 2003 to 
2013. The 2012/2013 overall number of in-patient dis¬ 
charges and deaths was 1,027,005, and that of day- 
patient was 516,127 in Hong Kong. Among all patients, 
33.4% were day-patients. 

To prevent nosocomial or healthcare-associated 
infections, especially airborne ones, hospital hygiene 
and infection control are necessary. The Centres for 
Disease Control and Prevention (CDC) provides guid¬ 
ance to help healthcare personnel to follow standard, 
contact, and airborne precautions when caring for hos¬ 
pitalized patients with known or suspected viral infec¬ 
tions. 2 Effective prevention measures, e.g. an airborne 
infection isolation room (AIIR), are especially crucial 
for control of acute respiratory infectious threats. 


Proper ventilation also plays a key role in infection 
control by minimizing airborne bacteria and viruses. 
Both the spread of severe acute respiratory syndrome 
coronavirus (SARS-CbK) during the largest nosoco¬ 
mial SARS outbreak in Hong Kong and the recent out¬ 
break of Middle East respiratory syndrome (MERS) in 
the South Korean hospitals revealed that airborne dis¬ 
ease transmission through inefficient hospital ward ven¬ 
tilation systems can lead to dire health consequences. 3-5 
For a balanced ventilation that delivers indoor air 
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quality and energy efficiency, general hospital wards 
have been designed to meet certain air change require¬ 
ments. 2,6-10 For instance, an air change rate (< ach ) ran¬ 
ging from 2 to 6h -1 is suggested to help decrease local 
mean age of air, 11 while 4h -1 is recommended for 
energy savings. ’ However, as new information 
becomes available, current air change requirements 
for hospital wards should be re-evaluated and updated. 

This study investigates the effectiveness of ventila¬ 
tion design strategies for general hospital wards in 
terms of virus removal capacity. A computational 
fluid dynamics (CFD) simulation of a typical general 
ward of Hong Kong hospitals was conducted and three 
respiratory viruses, namely MERS-CoF, SARS-CoF 
and H1N1 influenza virus, were chosen. The findings 
can be used by the hospital management to minimize 
cross-infection risk while reducing energy consumption 
of ventilation. 


Bioaerosol drag force 

The motions of spherical and non-spherical bioaerosols 
in a ventilated space can be calculated using equation 
(1) by integrating the force balance on the bioaerosols 
in terms of the drag force per unit particle mass per 
relative velocity F D (N s kg -1 m -1 ), where g is gravita¬ 
tional acceleration (m s -2 ), F x is the additional acceler¬ 
ation force per unit particle mass (N kg -1 ) if any, v b 
and v a are the velocities of the virus and air (m s -1 ), 
respectively, p a is the molecular viscosity of air 
(kgm -1 s -1 ), d b is the equivalent bioaerosol diameter 
(/xm), Re b is the Reynolds number for bioaerosols in 
an airflow field, p a is the air density (kg m -3 ) and p b is 

_O 1 A 

the virus density (=1100 kg vcT ). 


dvb 

dr 

F D : 


77 / \ - §(Pb Pa) r, 
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The drag coefficient C D for the bioaerosols is defined 
by equation (2) 


C D = : Re h < 1 (2) 

Re b 

Equations (1) and (2) are used in the CFD simula¬ 
tion to determine the motions of droplet nuclei under a 
Lagrangian scheme. 

The drag constant for the bioaerosols K D in equation 
(2) is given by equation (3) 


Validity of equation (3) was established for equivalent 
bioaerosol diameters d b = 0.69-6.9 /xm over a bioaerosol 
size range and further examined for Bacteriophage Phi 
X174 (ATCC 13706-B1) with d b = 0.054/xm. A combined 
experimental and simulation method described by Wong 
et al. 14 was performed. Before the experiment, agar plates 
were prepared in such a way that the entire preparation 
was covered by the host bacterium Escherichia coli 
(ATCC 13706). Plaque assays for bacteriophages were 
then conducted. The number of bacteriophages collected 
was measured by plate counting. 15 Figure 1 shows that 



Equivalent bioaerosol diameter db (jam) 


Figure 1 . Bioaerosol drag constant K D against equivalent 
bioaerosol diameter d h . 





Figure 2. Bioaerosol drag constants and absolute errors for 
Bacteriophage Phi XI74. 
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equation (3) can be extrapolated for predicting the disper¬ 
sion of smaller sized bioaerosols (down to d b = 0.054 /xm). 
When using K D = 0.001458 (Figure 2), the deviation 
between simulated and measured fractional counts 
inside a ventilated chamber was insignificant (p > 0.99, 
paired /-test). The corresponding fractional bias (FB) £j 
and normalized mean square error (NMSE) s 2 for 
K d — 0.001458 were 9.2% and —3%, respectively, as illu¬ 
strated in Figure 2. 16 

Using equation (4), where the projected image area, 
length and width are A b , / 7 and / 2 , respectively, the equiva¬ 
lent bioaerosol diameter d b can be determined from elec¬ 
tron micrographs. 17,18 Electron micrographs of samples 
of MERS-CoF, SARS-CoF, H1N1 influenza virus and 
Bacteriophage Phi X174 are shown in Figure 3. 


max (/ 1 , / 2 ) 

^ aspect — i 77 TT 

min (/i, / 2 ) 


CFD simulation 

A typical semi-enclosed six-bed general ward cubicle 
(7.5m (L)x6m (IF) x 2.7m (//)) employed for the 
CFD simulation is shown in Figure 4. As illustrated 
in Figure 5, there were four ceiling air inlets with 
mounted diffusers in the cubicle, and the corridor was 
the air-outlet. Four different air change rates (i.e. 
ach = 3, 6, 9 and 13h _1 ) were used, while the tempera¬ 
ture and relative humidity (RH) of the supply air were 
285 K and 80-95%, respectively. Six probable viral 
emission points, representing the positions of six 
patients lying on their respective beds (i.e. Man 1, 
Man 2,..., Man 6), were chosen as the source locations. 
The viruses investigated were MERS-CoF, SARS-CoF 
and H1N1 influenza virus. The metabolic rate of a 
reclining patient was assumed to be 0.8 MET (i.e. 
46.6 W m -2 ). 24 Flalf of the heat (23.3 W m~ 2 ) trans¬ 
ferred from patient skin surface by convection was 
assumed. Patient beds were set as rectangular box 
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Figure 3. Electron micrograph of reference viruses, (a) MERS-CoF 19 (b) SARS-CoF 20 (c) H1N1 influenza virus 21 and (d) 
Bacteriophage Phi XI74. 22 
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Figure 4. A six-bed general ward cubicle (dimensions in mm). 


for the worst scenarios that airflow under the beds was 
fully blocked by medical equipment and relation instal¬ 
lations. This arrangement may cause higher numbers of 
bioaerosol particle deposited on other patients. 

The simulation of airflow field and virus dispersion 
inside the ward cubicle was done through the use of 
CFD software FLUENT 14. An Eulerian-Lagrangian 
framework was used to solve the gas-solid two-phase 
flow problem, i.e. an Eulerian scheme for the prediction 
of a steady-state airflow field of the ward cubicle, fol¬ 
lowed by a Lagrangian approach for the determination 
of virus particle movements. 

In the Eulerian framework, a continuous phase of 
the induced airflow field was obtained with a second- 
order solution scheme. The renormalization group 


(RNG) k-s turbulence model was tested to be an appro¬ 
priate choice among the RANS turbulence models, and 
therefore, it was adopted to determine the air turbu¬ 
lence in the field. The model offered better accuracy 
and stability in cases of low Reynolds number and 
near-wall flows, and this model was found suitable for 

T A O ^7 

indoor airflow simulations also. ’ To couple with the 
pressure and velocity fields, the pressure implicit with 
splitting of operator (PISO) algorithm was employed, 
and the convection term was discretized using a second- 
order upwind scheme. 

To optimize the simulation quality and speed, three 
reference grid sizes namely fine (i.e. 2618 k), moderate 
(i.e. 1143 k) and coarse (i.e. 1071k) were constructed 
for determining a suitable mesh size. The fine gird 
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Figure 5. CFD configurations of a six-bed general ward cubicle. 


size was with a mesh skewness less than 0.25 to ensure 
an excellent cell quality. The moderate and coarse grid 
sizes were doubled and quadrupled by a second-order 
method for increasing the simulation speed. To assess 
the mesh quality using linear grid stretching, the 
asymptotic range of convergence c asymp applied 
was based on the grid convergence indexes (GCIs) 
with a relative error of computed average mass flow 
rate s rm s • 

If the value of c asymp in equation (5), where f s is the 
safety factor, r r is the refinement ratio and c Q is the 
theoretical order of convergence, is approximately 1, 
then the grid quality is fine. 


G Clf ine 

Casymp 

The GCI analysis results from the simulation were 
51 /o, 37 /o and 1.187 for GCI coarse , GCIj=i ne and c asvmp , 
respectively. 

In the CFD simulation, the transmission pathway 
of virus-laden respiratory droplets expelled by 
sneezing was predicted (i.e. d h initial virus-laden 
respiratory droplets diameter = 8.3 /xm, v^, = 50m s _1 , 
p b — 1100 kg m~ and n s , amount of virus particles 

1 A 9Q_O 1 

expelled by sneezing = 10,000 virus particles). ’ 
General patient wards are recommended at an RTI 



30-60%. 6 Within a short period of time (<0.1 s) after 
emission, the droplets would evaporate to droplet 
nuclei. 32 These nuclei were the dried-out residual of 
droplets possibly containing infectious pathogens. 
Only a small proportion (<10%) of the virus-laden 
droplets of a total number of 73,000-1,000,000 droplets 
expelled by a vigorous sneeze was assumed; neither 
aggregation with other particles nor cluster of virus 
would be formed from the dried-out nuclei at such 
low concentration. 34,35 Lipid-enveloped human corona- 
virus 299 E would remain alive by a half-life of 67 h at 
an RH of 50% and an air temperature of 20°C. 36 The 
survival rate of the three virus was expected to be 100% 
in simulation times less than 100 s. Other details for 
simulations and boundary conditions are summarized 
in Table 1. 

Table 2 exhibits the equivalent bioaerosol diameters 
d b for the virus droplet nuclei examined in this study. 
A one-way coupling was applied in the prediction to 
prevent the effect of particles on the continuum airflow. 
Each virus particle was tracked separately for its pos¬ 
ition and velocity by a discrete phase model (DPM). 
A previous study confirmed that the isotropic dis¬ 
crete random walk (DRW) model was effective and 
accurate in modelling bioaerosol dispersion and distri- 
bution due to turbulent fluctuations in the flow. For a 
coagulation effect of bioaerosol particles in this study, a 
very low volume fraction (<3000 cm - ) was kept in the 
ward cubicle to reduce collisions of the virus particles in 
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OO OQ 

a turbulent flow. ’ Both the amounts of virus par¬ 
ticles exhausted to the corridor n e and deposited on 
the ward room surfaces n d were counted. 40 A perfect 
sink boundary condition was applied to the ward room 
surfaces in order that the virus particles impinged onto 


Table 1. CFD simulations and boundary conditions. 


Computational domain 

7.5 m (L) x 6m (JV) x 2.7m (. H ), 
RNG k-s turbulence model, 
standard wall function. 

Mesh configuration 

1,143,766 cells 

Total supply airflow rate 

1.24 kg s- 1 (for ach = 3), 2.48 kg 
s- 1 (for ach = 6), 3.72kg s -1 (for 
ach = 9) and 5.37kg s -1 (for 
ach= 13), 285 K (air 
temperature) 

Each inlet airflow rate 

0.31kg s- 1 (for ach = 3), 0.62 kg 

(0.6 m x 0.6 m) 

s' 1 (for ach — 6), 0.93 kg s -1 (for 
ach = 9) and 1.34 kg s -1 (for 
ach = 13), 285 K (air 
temperature) 

Four diffuser 

Four-way spread-type, supply jets 

(0.6 m x 0.6 m) 

had an angle of 15° from ceiling, 
adiabatic and reflect boundary 
type. 

Corridor (6 m x 2.7 m) 

Pressure-outlet, 295 K (backflow 
temperature), adiabatic, escape 
boundary type. 

Walls, ceiling, floor and 

No slip wall boundary, 295 K 

beds 

(surface temperature), adia¬ 
batic, trap boundary type. 

Patients 

No slip wall boundary, 23.3 W 
m , trap boundary. 

Patient mouths 

Single-shot release with an upward 

(0.05 m x 0.05 m) 

exhalation velocity of v b = 50 m 
s" 1 , initial virus-laden respira¬ 
tory droplet diameter 
di— 8.3 pm, n s = 10,000 virus 
particles, density of bioaerosol 
particle p b = 1100 kg m -3 


RNG: renormalization group. 


the solid surfaces would be perfectly trapped with no 
reflection and desorption. 41 

Deposited ratio 

Direct contact with the viruses deposited on ward room 
surfaces and inhalation of the viruses suspended in the 
air are two potential transmission routes of virus-laden 
airborne particles expelled by sneezing 42 The deposited 
ratio r d is selected as a measure for evaluating the con¬ 
tact transmission; the movement of sneeze particles can 
be represented by equation (6), where r e is the exhausted 
ratio for particles exhausted to the corridor and r cl is the 
elapsed ratio for particles suspended in the air 

r e + r a + r d = 1 ( 6 ) 

Ratios r e , r a and r d are given by equation (7), where 
n s , n e , n a and n d are the amounts of virus particles 
expelled by sneezing, exhausted to the corridor, elapsed 



Figure 6. Description of the bioaerosol removal process in 
a general ward. 


Table 2. Virus information. 


Species 

ATCC 

Equivalent 
bioaerosol 
diameter d b (pm) a 

Aspect 

ratio 

(/aspect) 

Drag 

constant K D 

Evaporation 
time at 0% 

RH (s) 

Evaporation 
time at 50% 
RH (s) 

Evaporation 
time at 90% 
RH (s) 

MERS-CuF 

— 

0.167 ± 0.012 

1.27 

0.013945 

3.48 x 1(T 2 

9.55 x 10“ 2 

1.82 x 11T 1 

SARS-CuF 

— 

0.1375 ± 0.009 

1.052 

0.009453 

3.48 x 10“ 2 

9.55 x 1(T 2 

1.82 x 10“‘ 

H1N1 influenza virus 

— 

0.124 ±0.0001 

1.203 

0.0199 

3.48 x 10“ 2 

9.55 x 1(T 2 

1.82 x 10“‘ 

Bacteriophage Phi XI74 

13706-B1 

0.054 ±0.014 

1.012 

0.001458 

— 

— 

— 


MERS-CoF: Middle East respiratory syndrome coronavirus; SARS-CoF: severe acute respiratory syndrome coronavirus. 
Standard errors shown. 
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(a) Temperature (K) 
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Corridor 


Air velocity v a (m s ') 
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200801 
1 50e01 
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Air velocity v a (m s ') 
1 008+00 
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6.50801 
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500802 
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Figure 7. Simulation results of the ward with ach = 6 h . (a) Temperature distribution, (b) air velocity distribution and (c) 
flow pathlines from air supply inlets. 





Virus velocity 
v* (m s' 1 ) 

1 200^01 
1 90801 
1 00801 
1 70801 
1 60801 
1 50801 
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1 30801 
1 20801 




1 10801 
1 00801 
9 00802 
8 00802 
7 00802 
6 00802 
5 00802 
4 00802 
3 00802 
2 00802 
1 00802 
0 00e+00 


Figure 8. MERS-CoK pathways for six source locations with ach = 6 h \ (a) Man 1, (b) Man 2, (c) Man 3, (d) Man 4, 
(e) Man 5 and (f) Man 6. 
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(with respect to air suspension) and deposited on ward 
room surfaces, respectively 


Results and discussion 

Airflow field 


Figure 6 graphs the ratios r e , r cl and r d against time. 
The elapsed time x a (s) required for all virus particles to be 
exhausted or to deposit on ward room surfaces can be 
determined from the figure as represented by equation (8). 

r a (r a ) = 0 (8) 


Figure 7 shows the simulated temperature, air velocity 
distribution and the flow pathlines inside the ward with 
ach = 6h~ l . An increasing temperature gradient was 
found from the ceiling to the floor inside the ward. As 
shown in Figure 7(a), warmer air observed near the 
patients (due to body temperature) was diffusing out to 
the corridor. The average room temperature was about 
296 K. Stagnant air, with velocity less than 0.05 m s -1 , 
was found near the window side (i.e. Wall 2) and the 



Figure 9. Deposited and exhausted ratios of MERS-CuFfor six source locations, (a) Man 1, (b) Man 2, (c) Man 3, (d) Man 
4, (e) Man 5 and (f) Man 6. 
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lower floor level (i.e. below 0.3 m) as presented in 
Figure 7(b). The flow pathlines in Figure 7(c) confirm 
that there was insufficient ventilation to remove any 
virus particles in those areas. Generally, the airflow dir¬ 
ection was from the inner part of the ward to the 
corridor. 

Cross-infection from surface deposition 

Figure 8 indicates the potential risks of cross¬ 
infection with MERS-CoF through air pathways from 
different infected patients inside the ward. Patients stay¬ 
ing on the same side of an infected patient, especially 
the one located next to the corridor (i.e. Man 1 or 
Man 2), would have a higher chance of cross¬ 
infection. Two different virus pathway flows in the 
simulation due to the asymmetric diffuser locations 
are highlighted. Cases exhibited in Figure 8(c) and (e) 
show the virus moved along floor surface of the 
ward but in cases shown in Figure 8(d) and (f), virus 
would pass over nearby patients’ heads, then flew to 
the corridor. 

Figure 9 exhibits the values of deposited ratio r d and 
exhausted ratio r e of MERS-CoF on all surfaces inside 
the ward with ach = 3, 6, 9 and 13 h _1 . The results sug¬ 
gested that the deposition of virus particles was depend¬ 
ent on the location of an infected patient, and the virus 
particles would deposit mainly back on the source 
patient. If the infected patient was located near the cor¬ 
ridor (i.e. Man 1 or Man 2), the virus particles would 
likely be exhausted to the corridor. On the contrary, if 
the infected patient was located in the inner part of the 
ward (i.e. Man 5 or Man 6), the possibility that the 
virus particles would deposit on the wall surfaces or 
other patients was higher. 

Moreover, more virus particles would remain on the 
sneezing patient at lower air change rates, especially for 
ach = 6h~ l (r d >0.7 for all source locations), whereas 
more virus particles would deposit on the wall surfaces 
or be exhausted to the corridor at higher air change 
rates (i.e. ach = 9h _1 or 13h _1 ). Although the possibil¬ 
ity of cross-infection among patients and healthcare 
workers (HCWs) was lower at a higher ach (up to 
9h _1 ), HCWs should clean all patients and ward sur¬ 
faces regularly, regardless of ach. 

Cross-infection from inhalation 

Figure 10 graphs the MERS-CoF removal processes in 
the ward for the six source locations (i.e. Man 1, Man 
2,..., Man 6) and ach = 3, 6, 9 and 13h _1 . The results 
showed that virus particles from an infected patient 
located near the corridor were likely to be exhausted 
to the corridor. Meanwhile, lower r e (i.e. higher r d ) 
and longer elapsed time x a were associated with the 


MERS-CoF particles from the sneezing of an infected 
patient located in the inner part of the ward, and that 
indicated a higher risk of MERS-CoF infection 
through direct inhalation of particles or indirect inhal¬ 
ation of re-suspended particles. Besides, higher r e and 
shorter x a were found with increasing air change rates. 
If ach was increased from 3 to 13 h _1 , x a could be shor¬ 
tened by more than 30 s, and the risk of cross-infection 
from inhalation could be effectively lowered. Similar 
dispersion and deposition results were observed for 
SARS-CoF and H1N1 influenza virus particles and pic- 
torialized in Figures 11 and 12. The simulation per¬ 
formances demonstrated that virus particles with a 
relatively small equivalent bioaerosol diameter 
(d b ~ 0.1 /im) had similar particle dispersion and depos¬ 
ition characteristics in a general hospital ward. 

Figure 13 plots the average elapsed time x a against 
ach for MERS-CoF, SARS-CoF and H1N1 influenza 
virus. The results showed that x a could be significantly 
shortened by increasing the air change rate in the ward 
(R z >0.9). However, the threshold or optimal x a in the 
ward could not be determined as the infectious doses 
(ID50) of the three viruses varied in wide ranges (i.e. 
180, 1800 and 180 virus particles for MERS-CoF, 
SARS-CoF and H1N1 influenza virus, respectively, to 
cause a 50% infection). 43-45 Using the ASHRAE stand¬ 
ard for an AIIR (i.e. ach = 12h _1 ) as a safety measure, 6 
the corresponding elapsed time was about 34 s. When 
applying the ASHRAE (ach = 4h -1 ) and CIBSE 
{ach = 6 h -1 ) standards for a general patient room, the 
average values of x a were found to be 70s and 61s, 
respectively. 6,8 The elapsed time doubled when ach 
dropped from 12 to 4h _1 , and thus doubling the poten¬ 
tial inhalation risk. Based on the median value in 
accordance with both ASHRAE and CIBSE standards, 
the maximum ach in a general hospital ward should be 
9 h -1 (r a = 48 s) for the needs of maximizing energy effi¬ 
ciency and minimizing infection risk. Furthermore, it 
should be noted for other means to minimize the pos¬ 
sibility of cross-contamination in hospital wards, such 
as the installation of ultraviolet germicidal irradiation 
(UVGI) lamps for the destruction of viral nucleic 
acids. 46 


This study should provide a useful source of reference 
for the hospital management to mitigate the risk of 
infection with MERS or other airborne transmitted 
viruses through better ventilation design strategies. 
The results of this study demonstrated that the location 
of an infected patient would affect the infection risks to 
other occupants and HCWs inside the same hospital 
ward, and an increased air change rate in the ward 
could reduce the risk of infection from direct contact 


Downloaded from ibe.sagepub.com at DREXEL UNIV LIBRARIES on June 5, 2016 


10 


Indoor and Built Environment 0(0) 


e 

cd 


(N 

C 

c3 


c 


C 

G3 


c 

cd 


VO 

c 

c3 



0 20 40 60 80 100 


I 

0.8 

0.6 

0.4 

0.2 

0 






I 

0.8 

0.6 

0.4 

0.2 

0 



0 20 40 60 80 100 






0 20 40 60 80 100 


0 20 40 60 80 100 

ach = 6 h' 1 


ach = 3 h' 1 
x-axis: Time (s) 

y-axis: Deposited ratio rj, exhausted ratio r e 


I 

0.8 

0.6 

0.4 

0.2 

0 



0 20 40 60 80 100 





0 20 40 60 80 100 

ach — 9 h" 1 



20 40 60 80 100 


I 

0.8 

0.6 

0.4 

0.2 

0 





0 20 40 60 80 100 

ach =13 h" 1 


Figure 10. MERS-CoK removal processes for different source locations and air change rates. 
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Figure 11. SARS-CoK removal processes for different source locations and air change rates. 
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Figure 12. H1N1 influenza virus removal processes for different source locations and air change rates. 
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Figure 13. Average elapsed time z a with design standards. 

MERS: Middle East respiratory syndrome; SARS: severe acute 
respiratory syndrome; AMR: airborne infection isolation room. 


and inhalation. For a typical semi-enclosed six-bed gen¬ 
eral ward of Hong Kong hospitals, an air change rate 
of 9h _1 could effectively minimize the deposition and 
floating time of respiratory virus particles while max¬ 
imizing energy efficiency. In order to minimize the pos¬ 
sibility of cross-contamination in hospital wards, 
installation of UVGI lamps is also recommended. 
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